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Abstract. The Onsager formula for the free energy of the two dimensional Ising model with 
periodic boundary conditions is proved. 
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^ ■ 1 Introduction 

Onsager and Kaufman [||], [^, [Q proposed the formula for the partition function of the 

Qv^ ■ two dimensional Ising model. Another heuristic method for the calculation of the partition 

^ ■ function of the two dimensional Ising model is proposed by Kac and Ward [^]. Their idea 

'jrt . consists of the construction of the special matrix A whose determinant is connected with 

the partition function of the two dimentional Ising model. They considered simultaneously 

, _ two formulae: the determinant of the matrix A is proportional to the partition function of 

pH ■ the Ising model and it is proportional to the square of the partition function. For the proof 

O ■ of the first formula they used a topological statement. Sherman [^] constructed a counter - 

example for this statement. Hurst and Green [^ proposed to use for the calculation of the 

Ising model partition function not a determinant but a Pfaffian of some special matrix. This 

rS I method was improved in the papers |^, [§. The proof proposed by McCoy and Wu is 

c^ ' considered as the most mathematically correct. 

We consider a rectangular lattice on the plane formed by the points with integral Carte- 
sian coordinats x = j,y = k,0<j<N,0<k< M, and the corresponding horizontal and 
vertical edges in which the opposite sides of the entire rectangular are identified. In other 
words, we consider a rectangular lattice on a torus. McCoy and Wu proved the following 
formula for the partition function of the Ising model 

Z = (2cosh/5Eicosh/5E2)^^^l/2(-PfAi + Pfv42 + Pfv43 + PfA4), (1.1) 

where /3 is the inverse temperature and Ei{E2) is horizontal (vertical) interaction energy. 
McCoy and Wu could not calculate the Pfaffians Pf Aj. They calculated the determinants 
det Ai only. Since 

Pfi, = ±(detA,)^/^ (1.2) 
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it is necessary to calculate the sign in the formula ( |1.2| ). These signs are calculated in by 
using the heuristic method. 



In the paper |jTO[ a general formula of type ( |1 . 1|) is obtained for the lattices placed on 
an arbitrary orientable surface. In this paper we prove that for the torus case the general 
formula |T^ is the following formula 



Z = (-2cosh/5Ei cosh/?E2)*^^l/2(-Pf Ai + Pf A2 + Pf A3 + Pf A4), (1.3) 

We will calculate the Pfaffians PiAi and we will prove the Onsager formula for the free 
energy of the two dimensional Ising model. It is easy to guess our formulae for the Pfaffians 
by using the formulae for det Ai given in the book ^ and the equality 

(1 + zl){l + z^) - 2^i(l - zl) - 2^2(1 - zl) = 

{ZiZ2 + Z1 + Z2- if. (1.4) 

2 Homological Formula 

We denote by G{N, M) the graph described in the previuos section. Let the function a on 
the vertices of the graph G{N, M) takes the values in the multiplicative group Z2 = {1, —1}. 
The energy H{a) for the Ising model with zero magnetic field can be expressed in the form 

N M 

H{o) = - E E(^i^U k)a{j + l,k) + EMj, k)aij, k + 1)), (2.1) 

j=i k=i 

where the numbers 1 and A^ + 1 are identified and the numbers 1 and M + 1 are identified 
also. The partition function of the Ising model with zero magnetic field is defined as follows 

Z{N,M)=Y,exp{-PH{a)}, (2.2) 

(T 

where the summing runs over the group of the functions a on the vertices of the graph 
G{N, M) taking values in the group Z2. By performing the Fourier transformation |jll|, ([^], 
Chapter 5) we have 

2 
Z{N,M) = (2cosh/3Eicosh/5E2)*^^En(tanh/3Ei)'=^(«\ (2.3) 

The function ^ on the edges of the graph G{N, M) is called the cycle with the coefficients 
in the group Z2 if the product of all values of the function ^ on the edges connecting an 
arbitrary vertex of G{N,M) to other vertices is equal to 1. The number fci (0(^2(0) is the 
total number of all horizontal (vertical) edges on which the cycle ^ takes the value —1. 

By a dimer configuration on the graph G we mean a system of edges U = {{pk^Qk)} 
satisfying the following conditions: the edges of U are pairwise disjoint; all vertices of the 
graph G are covered by the edges of the dimer configuration U. 

Due to |T0| we will establish a one - to - one correspondence between the elements of the 



group of cycles with coefficiens in Z2 and the dimer configurations on another, expanded 
graph. Every vertex of the initial graph G{N, M) is connected with four other vertices. In 
the new graph we replace every vertex with two new vertices and add to the old four edges 



one new edge connecting two new vertices. We assume that new vertices are the end points 
of two old edges and one new edge. Thus every vertex is connected with exactly three other 
vertices. By applying this construction for every vertex of the graph G{N, M) we obtain a 
new graph tr{G{N,M)). Let a cycle ^ be given on edges of the graph G{N,M). Then it 
is given on the old edges of the graph tr{G{N,M)). By using the condition d^ = —1 we 
continue this function ^ as Z2 valued function ^ on the edges of the graph tr{G{N,M)). 
Namely for every vertex of the graph tr{G{N,M)) the product of the function ^ values 
on the edges connecting this vertex with three other vertices of the graph tr{G{N,M)) is 
equal to —1. The function ^ is a cycle and every vertex of the graph G{N, M) corresponds 
with two vertices of the graph tr{G{N, M)). Therefore the extension ^ on the edges of the 
graph tr{G{N, M)) is defined uniquely. Conversely the restriction of every function ^ on the 
edges of the graph tr{G{N, M)), satisfying the condition 9^ = — 1, to the edges of the graph 
G{N, M) defines a cycle on the edges of the graph G{N, M) with coefficients in the group 

Z2. 

We change every vertex of the graph tr{G{N, M)) by a triangle each vertex of which is 
an end point of one edge of the graph tr{G{N,M)) and an end point of the two edges of 
triangle. We obtain a new graph d{G{N,M)). The set of the functions ^ on the edges of 
the graph tr{G{N,M)), satisfying condition dC, = —1, is in one - to - one correspondence 
with the set of all dimer configurations on the graph cl{G{N,M)). Indeed, let us include 
into the dimer configuration U{^) all the edges on which the function ^ takes the value —1. 
Then three or one such edges come to every triangle. In the first case all vertices of the 
triangle are covered by the edges which do not intersect each other and some edges of the 
triangle does not belong to the dimer configuration U{^). In the second case one edge of the 
dimer configuraton f/(^) comes to one vertex of the triangle, and we include into the dimer 
configuration the opposite side of the triangle. Now again all vertices of the triangle are 
covered by the edges from the dimer configuration U{C,) which do not intersect each other. 

Let the dimer configuration U on the graph cl{G{N,M)) be given. Let us define the 
function C,{U) by taking it is equal to 1 on the edges of the graph tr{G{N, M)) which do not 
belong to the dimer configuration U and by taking it equal to —1 on the edges of the graph 
tr{G{N,M)) which belong to the dimer configuration U. If one side of a triangle belongs 
to the dimer configuration U, then one edge from the dimer configuration U comes to a 
triangle. If the sides of a triangle do not belong to the dimer configuration U, then three 
edges from the dimer configuration U come to a triangle. In both cases the function ^{U) 
satisfies the condition d^ = —1. 

One - to - one correspondence allows us to rewrite the partition function ( p.3|) in another 
form. Every dimer configuration U on the graph cl{G{N, M)) corresponds with a cycle ^ 
on the edges of the graph G{N, M) with coefficients in the group Z2: (p, g) G f/ fl G{N, M) 
if and only if C,{{p,q)) = —1. Therefore the relation ( |2.3| ) can be rewritten in the following 

form 

2 

Z{N,M) = (2cosh/?Eicosh/5E2)*^^5]n(tanh/?Ei)#(^"^(^'*'^»% (2.4) 

U i=l 

where summation runs over all dimer configurations U on the graph cl{G{N,M)). ^{U fl 
G{N,M))i {^{U n G{N,M))2) denotes the total number of horizontal (vertical) common 
edges in the dimer configuration U and in the graph G{N, M). 

With each vertex of the graph ^(A^, M) there corresponds a cluster of the graph 
d{G{N, M)): two triagles connected by one edge. We place one triangle under another such 
that their bases are horizontal. Let vertex opposite to the base of the upper triangle be set 



below the base and the vertex opposite to the base of the lower triangle be set above the 
base. The vertical edge connects these two vertices. The right vertex of the upper triangle 
is denoted by 1. The left vertex of the lower triangle is denoted by 2. The left vertex 
of the upper triangle is denoted by 3. The right vertex of the lower triangle is denoted 
by 4. The upper vertex of the lower triangle is denoted by 5. The lower vertex of the 
upper triangle is denoted by 6. With all six vertices there correspond two numbers {j,k), 
^'^j^N,l<k<M corresponding to the old vertex of the graph G{N, M). The vertices 
(l,j, /c) and (6,j, /c); (l,j, /c) and (3,j, /c); (3,j, fc) and (6,j, /c); (5,j, /c) and (6,j, /c); (2,j, /c) 
and (4,j, fc); {2,j,k) and (5,j, fc); (4,j, fc) and (5,j, /c) are connected by the edges of one 
cluster. The vertices (1, j, k) and (2, j + 1, k) are connected by the horizontal edges of the 
graph G{N, M). We consider the numbers 1 and A^ + 1 identical. The vertices (3, j, k) and 
(4, j, k + 1) are connected by the vertical edges of the graph G{N, M). The numbers 1 and 
M + 1 are considered identical. 

Let an orientation of every edge of the graph cl{G{N, M)) be given. If an edge is defined 
by its end points {p,q), then an orientation is given by a function (l){p,q) taking only two 
values 0, 1 and satisfying the following condition 

0(p, g) + 0(g,p) = 1 mod2. (2.5) 

Let us define the function z{p, q) = z{q,p) on the oriented edges (p, q): z{p, g) = 1 if the 
vertices p,q belong to the same cluster and are connected by an edge; z{p,q) = tanh/3i5i 
{z{p, q) = tanh/5i?2) if the vertices p, q belong to the neighboring clusters and are connected 
by a horizontal (vertical) edge; z{p, g) = if the vertices p, q are not connected by an edge. 

Let the linear ordering of vertices of the graph cl{G{N, M)) be given 

n{i,j, k) = t + 6{k-l) + 6M(j - 1) (2.6) 



where l<«<6,l<j<A^, 1<A;< M. The function ( p.6|) establishes the correspondence 
between the vertices of the graph cl{G{N, M)) and the numbers 1, ..., QMN. 
Let us define 6MN x QMN - matrix 

A{(j))n(p)n{q) = exp{m0(p, q)}z{p, q) (2.7) 

where the numbers p = {i,j, k) are indices for the vertices of the graph cl{G{N, M)). Due 
to the relation (|2.5|) the matrix (|2.7| ) is antisymmetric. Hence we can define its Pfaffian 



Pf A(0) = ((3MiV)!)-^2-3^-^^ Y. (-1)"^"^^(0).(1K(2) ■ ■ ■ ^(0).{6MiV-l).{6MiV) (2.8) 

where vr is an arbitrary permutation of the numbers 1, ..., QMN and (t{'k) is its parity. 
The definition ( |2.8| ) can be expressed as 

pf ^(0) = Ei-ir^""^ n ^(</')m (2-9) 

U (i,j)<^U 

where U is an ordered subdivision of the numbers 1, ..., QMN into pairs 
{(H,J2),---,(»i3M7v,j3AfAr)}, 1 = H < ■■■ < Hmn < QMN, ip < jp, p = 1,...,3MN. 
The number (t{U) = 0, 1 respectively with the parity of the permutation mapping the 
numbers {ii,ji,...,i3MN,J3MN) into the numbers (1, ..., 6MA^). Due to the relation ( |2.71 ) 



and the definition of the function z{p, q) the summation in ( p.9| ) runs only such subdivi- 
tions U which correspond with the dimer configurations on the graph cl{G{N,M)). Let 
U = {{pi,qi), ■■.,{p3mn,(13Mn)} be a dimer configuration on the graph d{G{N,M)). By 
using an orientation function satisfying the condition ( |2.5| ) we define a function on the 
dimer configurations on the graph d{G{N,M)) 

3MN 

r[0](f/) = a{n{pi),n{qi), ..., n(p3M7v), ?^(g3M7v)) + J2 <^(Pi'?i) "^0^2 (2.10) 

where a{n{pi), n{qi), ...,n{p3MN), n^q^MN)) = 0, 1 respectively with the parity of the permu- 
tation mapping the numbers (n(pi), n{qi), ...,n{p3MN),n{q3MN)) into the numbers (1, ..., 6MN). 
The function 0(p, q) satisfies the condition (|2.5|). Hence the the right hand side of the equal- 
ity ( |2.1Q| ) is independent of an ordering of the vertices in the edges {pj,qj) although every 



summand in ( p.lO|) depends on this ordering. The right hand side of the equality (|2.10|) is 
also independent of an ordering of the edges {pj, qj). 

The substitution of the definitions (|2.7] ) and (|2.1CI| ) into the right hand side of the equality 



(^ yields 

Pf A(0) = ^(-l)^[^](^);z([/) (2.11) 

u 

where 

ziU)= n <P^l) (2-12) 



and the summation in ( p.ll| ) runs over all dimer configurations on the graph cl{G{N, M)). 



The expression ( |2.11|) differs from the sum ( p.4|) by the multiplier (— l)'^M(t^) only. For 
any dimer configuration Uq on the graph cl{G{N, M)) we have 

Pf A(0) = (-l)-['^K^o) ^(_l)rM(f/o)+rM(t/)^(^)_ (2.13) 

u 

Let us study the function r[0](f/i) +r[0](f/2) of a pair of the dimer configurations Ui, U2 
on the graph cl{G{N,M)). Let us consider the symmetric difference 

f/iAf/2=(f/iUf/2)\(f/inf/2). (2.14) 

The dimer configuration definition implies that this set consists of the finite number of closed 
broken lines which do not intersect each other, that is, we can write 

f/iAf/2 = UU^i. (2.15) 

Every broken line 7^ does not intersect itself. The edges from Ui and U2 are included in each 
broken line alternatively. Due to Theorem 2 from |iy] for any dimer configurations f/i, f/2 
on the graph d{G{N, M)) 

rmU^)+T[ct>\{U2) = s + j2(t>{l^) mod2 (2.16) 

where the broken lines 7^ are given by the relation ( [2.15| ) and 

0(7i)= E 0(p,g)mod2. (2.17) 



This sum does not depend on the orietation chosen on the broken hne 7^ since the total 
number of edges in the broken hne 7^ is even. 

Let us define on the graph cl{G{N,M)) special orientation function 

<l){{l,j,k),{6,j,k)) = <j){{6,j,k),{3,j,k)) = 
<P{{3,j,k),{l,j,k)) = 0mod2 
^{{2,j,k),{5,j,k)) = 4>{{5,j,k),{4,j,k)) = 
<P{{A,j,k),{2,j,k)) = 0mod2 
0((5,j,A;),(6,j,A;)) = Omod2 
cl>{{l,j,k),{2,j + l,k)) = 0mod2 
cf>{{3,j,k),{4,j,k+l)) = 0mod2. (2.18) 



In the equalities (|2.18|) the numbers 1 and A^ + 1, 1 and M + 1 are identified. The graph 



cl{G{N^ M)) consists of the triangles and 12 - angles. The relations ( ^.1^ ) have the property: 
if a broken line 7 is a boundary of a triangle or 12 - angle and the counter - clockwise 
orientation is chosen, then 

0(7) = 1 mod 2. (2.19) 

Theorem 2.1 Let on the graph cl{G{N,M)) the orientation function (|2.18|) be given. Let 
domain on the torus be bounded by the broken lines 71, ...,7^ and these broken lines coincide 
with the symmetric difference of two dimer configurations on the graph cl{G{N, M)). Then 

s 

s + 5l0(7i) = 0mod2. (2.20) 

Proof. In view of the condition of the theorem the domain on the torus is bounded by 
broken lines 71, ...,7^ which do not intersect each other. We stretch the auxiliary disk di on 
each broken line 7^. In general, we obtain the new closed orientable surface P. The domain 
consists of the faces a (triangles and 12 - angles). We choose the orientation which induces 
the counter- clockwise orientation on the broken lines da. Let us choose an orientation of 
each auxiliary disk di coherently with the orientation of the boundary of the domain. 

Since the set of closed broken lines 71, ..., 7^ coicide with the symmetric difference of two 
dimer configurations, the domain contains an even number of vertices, the number of edges 
and the number of vertices for each broken line 7j are also even. Therefore the number of 
vertices v of the constructed cell complex for an orietable closed surface P is even. Let / 
and e denote the total numbers of faces and edges of the orientable surface P. By using the 
skew symmetry condition ( p.5|) and by summing over all faces of the surface P including the 
auxiliary disks d^ we have 

^(1 + 0(9(t)) = / + e = t; + x = 0mod2, (2.21) 

since Euler characteristic x = 2(1 — g) of the orientable surface P is even. 

If the orientation function on the graph cl{G{N,M)) is given by the relations ( p.l8| ). 



then for any face of the surface P excepting the auxiliary disks the relation ( p.l9|) holds. 
Hence the left hand sides of the equalities ( p^.20D and ( |2.21| ) coincide. The theorem is proved. 



Theorem 2.2 Let the dimer configuration Uq on the graph cl{G{N,M)) correspond to the 
cycle ^0 with the coefficients in the multiplicative group Z2 which takes the value 1 on all 

6 



edges of the graph G{N, M). Let the function T[(f)]{U) be given by the relation ( p.lO| ) for the 
orietation function defined by the equalities ( [2.181 ). Then 

r[0](?7o) = MN mod 2. (2.22) 

Proof. The dimer configuration Uq does not contain the edges connecting different clusters. 
In every cluster (j, fc) the dimer configuration Uq contains three edges: ((1, j, fc), (3, j, A;)), 



((2, j, k), {A J, k)) and ((5, j, k), (6, j, k)). In according to (U) 



a(n(l,l,l),n(3,l,l),n(2,l,l),n(4,l,l),n(5,l,l),n(6,l,l),..., 
ri(l, N, M), n{3, N, M), n(2, A^, M), n(4, A^, M), n(5, A^, M), n(6, N,M)) = 

MN mod2 (2.23) 

since there is one disorder: n(3, j, fc), n{2,j, k) in any of MN clusters. 
The definitions (p.l8|) and the relation ( p.5|) imply 



N M 

^J2Ml,j,k),{3,j,k))+<j){{2,j,k),{AJ,k)) + <t){{5J,k),{6,j,k))]=0mod2. (2.24) 
j=i k=i 



The substitution of the equalities (|2.23|) and (|2.24| ) into the equality (|2.1(]| ) for the dimer 



configuration Uq yields the equality ( |2.22 ). The theorem is proved. 

The number (|2.22|) is the total number of the vertices of the graph G{N, M) or it is the 
total number of the clusters of the graph tr{G{N, M)). The formulae ( p.. I]) and ( |1.3|) differ 
each other by the multiplier (— l)^[<^](^o)_ 

We denote by Z'^'^'^ the group of modulo 2 residuals. The modulo 2 residuals are multiplied 
each other and the group Z'^'^'^ is a field. The graph cl{G{N, M)) divides the torus into the 
faces: the triangles and 12 - angles. The cell complex P{cl{G{N,M))) = P{G) is called 
the set consisting of the cells (vertices, edges, faces). To every cell sf there corresponds the 
natural number p (dimension). For the vertices p = 0, for the edges p = 1 and for the faces 
p = 2. To every pair of the cells sf , Sj~ there corresponds the number (sf : s^" ) G Z'^'^'^ 
(incidence number). If the cell s^~ is included into the boundary of the cell sf, then the 
incidence number (sf : s^~ ) = 1. Otherwise the incidence number (sf : s^~ ) = 0. For any 
pair of the cells sf , s^ the incidence numbers satisfy the condition 

Y^{sl.sl)isl:s'^) = 0mod2. (2.25) 

m 

Indeed, if the vertex s^ is not contained in the boundary of the triangle or 12 - angle sf, then 
the condition ( p.25|) is fulfilled. If the vertex s° is included into the boundary of the face sj, 



then it is included into the boundaries of three edges s^ two of which are included into the 
boundary of the face s"^. The condition ( p.25|) is fulfilled again. 



A cochain c^ of the complex P{G) with the coefficients in the group Z'^'^'^ is a function 
on the p dimensional cells taking values into the group Z'^'^'^. Usually the cell orientation is 
considered and the cochains are the antisymmetric functions: c^(— s^) = — c^(s^). However, 
— 1 = 1 mod 2 and we can neglect the cell orientation for the coefficients in the group Z'^'^'^. 
The cochains form an Abelian group 

(cP + c'P)(sf) = cP{sf) + c'P(sf) mod 2. (2.26) 



It is denoted by Cp{P{G),Z^'^'^). The mapping 

dc^sr') = E(4 ■■ ^r')c(4) ^od2 (2.27) 



J 



defines the homomorphism of the group Cp{P{G), Z'if^) into the group C^ ^{P{G), TI^'^). It 
is called the boundary operator. The mapping 

9V(sr') = YX^'t^ ■ 4)c(4) mod 2 (2.28) 

j 

defines the homomorphism of the group C^(P(G'), 23'^'^) into the group C^+^(P(G'), Zj''''). 
It is called the cobounbary operator. The condition ( p.25|) implies dd = 0, d*d* = 0. 
A kernel Zp{P{G),Zf'^) of the homomorphism d : GP{P{G),Zf'^) -^ CP-\P{G),Zf'^) is 
called a group of cycles of complex P{G) with the coefficients in the group Z'^'^'^. The image 
Bp{P{G),Zf'^) of the homomorphism d : GP+\P{G),Zf'^) -^ GP{P{G),Zf'^) is called 
a group of boundaries of the complex P{G) with the coefficients in the group Z'^'^'^. Since 
dd = 0, a group Bp{P{G), Z'^'^) is a subgroup of the group Zp{P{G), Z'^'^). A quotient group 
Zp{P{G),Zf'^)/Bplp{Gizf'^) is called a homology group Hp{P{G),Zf'^) of the complex 
P{G) with the coefficients in the group Z'^'^'^. Similarly, for the coboundary operator d* a 
group of cocycles Z'p{P{G\Z^'^\ a group of coboundaries i?*'(P(G), Zj*^*^) and a group of 
cohomologies Hp{P{G), Z^'^) are defined. 

It is possible to introduce the bilinear form on Gp{P{G), Z'^'^^ 

(r,^^) = E/'(^?M^Dmod2. (2.29) 

i 

The definitions f^J^ and (g]^ imply 

(r,5V"') = W^^^-')mod2 (2.30) 

(/^ 9(7^+1) = (97^(7^+l) mod 2. (2.31) 

We identify every set of cells {s^} with a cochain taking value 1 on these cells and which 



takes value on other cells. For example, the symmetric difference ( p.l5|) of two dimer 
configurations IJ\^ U2 on the graph cl{G{N,M)) is identified with the cochain which takes 
the value 1 on the edges belonging to the broken lines 7^ and which takes the value on the 
other edges of the graph d{G{N, M)). 

In according to the relation (|2.13| ) it is necessary to study the function 

tpiUAUo) = r[(j)]{U) + T[(f)]{Uo) mod 2 (2.32) 

where the function r[0]([/) is defined by the relation ( ^.101 ) by using the orientation function 
( p.l8|) . For the dimer configuration Uq we take one corresponding to the cycle ^0 with the 



coefficients in the group Z2 which takes value 1 on all edges of the graph G{N,M). (We 
consider on the graph G{N, M) the cycles with the coefficients in the multiplicative group 
Z2. We hope that making use of the cycles with the coefficients in the different groups Z2 
and Zg'^'^ on the different graphs does not put us to cofusion.) Due to relations ( |2.15| ), ( |2.16| ) 



UAUo = Y.l^ (2-33) 



i=l 



V'([/A[/o) =s + ^7imod2. (2.34) 

Theorem 2.3 The function i/j(U AUq) , given by the relations ( |2.32| ) - ( p.34| ), depends on the 
equivalence class of the cycle UAUo in the homology group Hi{P{G),Z2'^'^) only. 
Proof. Due to relation ( p.33| ) the symmetric difference UAUq is the set of the broken 
closed lines. Hence it is identified with the cycle from the group Zi(P(G'), Zg'^'^). Thus 
the symmetric difference UAUq generates some equivalence class in the homology group 
ifi(P(G'), Zg'^'^). Let for the dimer configurations f/i and f/2 on the graph cl{G{N, M)) the 
symmetric differences UiAUq and U^AUq generate the same equivalence class in the ho- 
mology group. Let us prove ipiUiAUo) = iP{U2AUq). Indeed in this case the symmetric 
difference {UiAUq) A{U2AUq) is a boundary of some domain on the torus. Due to Lemma 
1 from 13 

(f/iAf/o)A(f/2Af/o) = U1AU2. (2.35) 

Therefore the symmetric difference U1AU2 is a boundary of some domain on the torus. Now 
the definition ( |2.32| ) implies 

^(U^AUo) + ijiU2AUo) = r[0](f/i) +r[0](f/2) mod 2. (2.36) 

It follows from the equalities ( |2.15| ), ( p^.l6D and equality ( |2.2(J| ) from Theorem 2.1 that the 

right hand side of the equality ( p.36| ) equals 0. Thus the values ipiUiAUo) and 'iIj(U2AUq) 

coincide. The theorem is proved. 

Theorem 2.4 The symmetric difference UAUq of the dimer configurations U, Uq on the 

graph cl{G{N, M)) generates all elements of the homology group Hi{P{G)), Z'^'^'^). 

Proof. Let the dimer configuration Ua consist of the edges ((3, j, 1), (6, j, 1)), ((4, j, 1), (5, j, 1)), 

J = 1,...,N; ((l,j,l),(2,j + l,l)), J = l,...,iV-l; ((1, AT, 1), (2, 1,1)); {{l,j,k),{3,j,k)), 

((2,j,A;),(4,j,A;)), ((5, j, A;), (6, j, fc)), l<3<N,2<k<M. Let us remind that the 

dimer configuration Uq corresponding to the cycle ^0 with the coefficients in the group Z2 

which takes value 1 on all edges of the graph G{N, M) consists of the following edges of 

the graph c/(G(iV,M)): ((1, j, fc), (3, j, A;)), ((2, j, A:), (4, j, fc)), ((5, j, fc), (6, j, fc)), 1 < j < iV, 

1 < fc < M. Thus the symmetric difference UaAUo is the closed non - intersecting itself non 
- homological to zero broken line which goes around the torus horizontally. 

Let the dimer configuration Uh consist of the edges: ((1, 1, k), (6, 1, k)), ((2, 1, k), (5, 1, k)), 
k = l,...,M; {{3,l,k),{A,l,k + l)),k = l,...,M-l; ((3, 1, M), (4, 1, 1)); ((1, j, A;), (3, j, A;)), 
((2, J, A;), (4, j, A;)), ((5, j, A;), (6, j, A;)), 2 < j < A^, 1 < A; < M. Thus the symmetric difference 
UbAUo is the closed non - intersecting itself non - homological to zero broken line which goes 
around the torus vertically. 

Let the dimer configuration Ua+b consist of the edges: ((5, 1, 1), (6, 1, 1)); ((3, j, 1), (6, j, 1)), 
((4, J, 1), (5, J, 1)), 2<j<N; ((1, 1, A;), (6, 1, A;)), ((2, 1, k), (5, 1, A;)), 2 < k < M; 
((l,j,l),(2,j + l,l)),j = l,...,iV-l;((l,iV,l),(2,l,l));((3,l,A;)(4,l,A;+l)),A; = l,...,M- 
1; ((3,1,M),(4,1,1)); ((1, j, A;), (3, j, A;)), ((2, j, A;)(4, j, A;)), ((5, j, A;), (6, j, A;)), 2 < j < iV, 

2 < k < M. Thus the symmetric difference Ua+bAUo is the closed non - intersecting itself 
non - homological to zero broken line which goes around the torus horizontally and vertically. 

At last UqAUq = 0. The equivalence classes of the symmetric differences UqAUq, UaAUo, 
UbAUo, Ua+bAUo as the closed curves on the torus represent all elements of the homology 
group Hi{P{G), Z'^'^'^). The theorem is proved. 

Let us define the cochain a* G G^{P{G), Z'^'^^ in the following way: 



a*{{{l, N, k), [2,1, k))) = 1, 1 < k < M. The cochain a* equals on all other edges of the 
graph cl{G{N, M)). By making use of the definition ( |2.28D it is easy to verify d*a* = i. e. 
a* G Z^(P{G), TJ^'^). Due to definition (|2.29|) the linear function (a*, /) is given on the group 
Zi{P{G),Zf'^). The equality (g^^) and the equality d*a* = imply that this function is 
equal to zero on the group Bi{P{G),Z2'^'^). Hence the linear function {a*,f) defines the 
function on the homology group Hi{P{G), Z'^'^'^). We denote by UaAUo the cochain which 
takes the value 1 on all edges of the broken line UaAUo. It is equal to zero on all other 
edges of the graph d{G{N, M)). The cochains UbAUo, Ua+b^Uo and UqAUo are defined 
analogously. The definition ( p^.29[ ) implies 

{a*, UaAUo) = 1 mod 2, {a*, UbAUo) = mod 2, 
{a*, Ua+bAUo) = 1 mod 2, (a*, UqAUq) = mod 2. (2.37) 

We define the cochain b* E C^{P{G), Zf'^) in the following way: 6*(((3, j, M), (4, j, 1))) = 
1, 1 < j < A^. The cochain 6* equals on all other edges of the graph cl{G{N,M)). By 
making use of the definition ( |2.28| ) it is easy to verify d*b* = i. e. b* e Z^{P{G),Z'!f'^). 
Due to definition ( p.29|) the linear function (6*, /) is given on the group Zi(P{G), Z^'^). The 
equality ( |2.31| ) and the equality d*b* = imply that this function is equal to zero on the 
group Bi[P{G), Z'^'^\ Hence the linear function (6*, /) defines the function on the homology 
group Ex{P{G),Zf'^). The definition i ^^^ implies 

(r, UaAUo) = mod 2, {b*, UbAUo) = 1 mod 2, 
{b*, Ua+bAUo) = 1 mod 2, {b*, UqAUq) = mod 2. (2.38) 



Let the orientation function satisfy the skew symmetry condition (p.5|) . Let us define 
three new orientation functions on the directed edges 

(0 + a*)((p,g))=0(p,g) + a*(p,g)mod2, (2.39) 

(0 + 6*)((p,g))=0(p,g) + 6*(p,g)mod2, (2.40) 

{^ + a* + b*){{p,q)) =(f){p,q) + a*{p,q) + b*{p,q) mod2. (2.41) 

Since a*{q,p) = a*{p,q), b*{q,p) = b*{p,q), the orientation functions + a*, (p + b*, (f) + a* + b* 



satisfy the skew symmetry condition ( p.5[ ). 

Theorem 2.5 Let the orientation function ( p.l8| ) be given on the graph cl{G{N,M)). The 



relation i\2.1\ ) defines four antisymmetric QM N xQMN - matrices: A((f)), A{(f)+a*), A{(f)+b*), 
A{(f) + a* + b*). Then the following relation for the partition function (p.2|) holds 



Z{N,M) = {-2 cosh f3Ei cosh pE2)^^^l/2[-PiA{4>) + Pi A{4> + a*) + 

FiA{(j) + b*)+FiA{(j) + a* + b*)]. (2.42) 

Proof. The relation ( p.4|) may be rewritten as 

Z{N, M) = (2 cosh (3Ei cosh /JEa)^^^ ^ z{U) (2.43) 



where the monomial z{U) is given by the equality (p.l2|) and the summation in (|2.43| ) runs 
over all dimer configurations on the graph cl{G{N, M)). 
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The equalities ( |2.1CI| ), (|2.18| ) and ( p.32|) define the function i(j{UAUo) of the symmetric 
difference UAUo of the dimer configurations U, Uq on the graph cl{G{N, M)). The dimer con- 
figuration Uq corresponds with the cycle ^o with the coefficients in the group Z2 which takes 
the value 1 on all edges of the graph G{N, M). Due to Theorem 2.3 the function il){U /\Uq) 
depends on the equivalence class of the cycle UAUq in the homology group Hi{P{G), TUf'^^ 
only. Hence it is sufficient to calculate the function iP{UAUq) on the following symmetric 
differences: UaAUo, UbAUo, Ua+bAUo and UqAUq. It is easy to verify 



V>(f/,Af/o) = ^(f/;,Af/o) = ij{Ua+bAUo) = 1 mod2,^(t/oAf/o 



0mod2. 



^2.44) 



By making use of the equalities ( |2.37| ), ( p.38| ) and ( p.44| ) it is possible to verify the following 
relation 

for four symmetric differences: UaAUo, UbAUo, Ua+bAUo, UoAUo- All terms of the equal- 
ity ( p.45| ) depends on the equivalence class of the cycle UAUo in the homology group 
Hi{P{G),7a2'^'^) only. Hence if the equality ( 2.45 ) holds for the symmetric differences: 
UaAUo, UbAUo, Ua+bAUo, UoAUo, it holds for the symmetric difference UAUo where U 
is an arbitrary dimer configuration on the graph cl{G{N, M)). 

For any dimer configuration U on the graph cl{G{N, M)) we introduce the cochain U G 
C^{P{G), 7/^^'^) which equals 1 on all edges from the dimer configuration U and which equals 
on all other edges of the graph cl{G{N, M)). Due to definition the cochains a* and h* are 
equal to zero on all edges from the dimer configuration f/g. Hence 



{a*, UAUo) = {a*, U) mod 2, {h* , UAUo) 



U) mod 2. 



(2.46) 



Let us multiply by (— l)'^(^^^o) the equality ( p.45| ). Then it follows from the equalities (p.32| ), 
( p.46|) and from Theorem 2.2 that 



\MN 



l/2[-( 



vmiu) 



+ (-1) 



T[4,](U)+{a*,U) 



+ 



(-1) 



T[4,](u)+{h* ,u) 



+ (-1) 



T[4>]{U)+{a\U) + {h*,U)l. 



(2.47) 



Now substituting the relation ( |2.47| ) into the equality ( f^.43D and making use of the equahties 
( p.7|) , (|2.10|) , (p.ll|) we obtain the equality (|2.42|) . The theorem is proved. 

For the calculation of the Pfaffians it is necessary to define the matrices A(0), A{(f) + a*), 
A{(p + h*), A{(p + a* + h*) exphcitly. The relations (^ - (^ and (|2iq ) imply 



A{4>)n{l,nM)n{2,J2M) — -^(0)n{2,J2,fe2)n(l,ii,fci) — 

tanh/?Ei((5ji+ij2 + 5j^N5j2i)hik2 

^(0)n{l,ji,A:i)n{3,J2,fc2) = ■~^(0)n(3,i2,fe2)n(l,ii,A:i) = -^juJ2^kik2 
^(0)n(lji,fci)n(6,i2,fc2) = ~^{(P)n{6,J2,k2)n(lJiM) = ^jld2^kik2 

A{(p)n{2,jiM)n(4:,J2,k2) = "~^(0)n(4,i2,fc2)n(2ji,A:i) = -^ji,J2^kik2 
M4>)ni2,jiM)ni5,J2,k2) = -^(0)n{5,J2,fc2)n{2,ji,fci) = ^jij2^kik2 

^{4>)n(3,ji,ki)n(4:,J2,k2) = ■~^(0)n{4,J2,fc2)n(3,ii,fci) = 

tanh/3ii^2'^jij2('^fci+i,A:2 + (^feiMi^tei) 



(2.48) 

(2.49) 
(2.50) 
(2.51) 
(2.52) 

(2.53) 
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^(</')n(3,ji,fci)n(6j2,fc2) = ~M't^)n{e,J2,k2)n{3,ji,ki) = ~^jl,J2^kik2 (2.54) 

^{(P)n{4,,ji,ki)n{5,j2,k2) = ~^i4>)n(5,J2,k2)n(4:,ji,ki) = ~^ji,J2^kik2 (2.55) 

^(0)n(5,ii,fci)n(6,i2,fc2) = ^^(0)n(6,J2,A:2)rt(5Ji,fci) = ^ji,J2^kik2 (2.56) 

All other elements of the matrix A{(j)) are equal to zero. 

In view of the definition the cochain a* equals 1 on the edges ((1, A^, k), (2, 1, k)), 1 < k < 
M and it is equal to zero on all other edges. The equalities ( |2.5| ) - ( |2.7| ) and ( p.l8| ) imply 



A(0 + a*)n{lj-i±j^)n{2,J2,k2) ~ ~^i^ + (^*)n{2,J2,k2)'n.{l,ji,ki) ~ 

tanh/5-E'i((5j^+ij2 — SjiNSj2i)^kik2 



(2.57) 



All the other elements of the matrices A{(j)) and A{(j)+a*) coincide. These non - zero elements 
are given by the relations ( p.49[ ) - ( |2.56| ). 

Due to definition the cochain b* equals 1 on the edges ((3, j, M), (4, j, 1)), 1 < j < A^, 
and it is equal to zero on all other edges. The equalities ( ^.51 ) - ( |2.7| ) and ( |2.18D imply 



A{(f) + b*)n(3j^^ki)n{4,J2M) " ~M^ + b*)n{4:,J2M)n{3,jiM) " 

tanh l3E26j-^jr,{6k^+i^k2 — hiMh2i] 



(2.5^ 



All the other elements of the matrices ^4(0) and A{(j)+b*) coincide. These non - zero elements 
are given by the relations (g^ - (|23^ ), (|23^ ) - (p36| ). 

In view of the relations ( p.5|) - ( p.7|) and (|2.18|) and the definitions of the cochains a*, b* 
we have 

A{(f) + a* + b*)n(l,j-^,ki)n{2,J2,k2) = ~M4> + 0.* + b*)n(2,J2M)n(l,hM) = 

tanh/3Ei(5j^+ij2 - 5j^N5j2i)^k^k2 (2-59) 



A{(t) + a* + b*)n{3,n,ki)n{'i,J2,k2) = ~M(P + a* + V')n{i,J2,k2)n{3,hM) = 

tanh /3E26j^J2{6ki^i,k2 — hiMh2i) 



(2.60) 



All the other elements of the matrices A{(f)) and A{(j) + a* + b*) coincide. These non - zero 
elements are given by the relations {^^ - ([23^ ), (P3^ ) - ( P3B|) . 

In the next section we calculate Pfaffians of the matrices (|2.48| ) - ( p.60| ). 



3 Pfaffians 

For the natural numbers — A^ < J < A^ the following relation is valid 

N 

i'=i 



The relations (|2^) - (p36D and {^ imply 



N M 



^W)n{ii,ji,ki)n{i2,J2,k2) 2-^ 2-^ / . 

H'*2~-'- ■?! ^2"-'- k-^,k2=l 

-S(0)„(i;jJ,fc;)n(nji,fci)^(0)n(i;j;,fci)n(i^j^,fc^)5(0)n(i^,j^,fc^)n(i2,J2,fc2) 



(3.1) 



(3.2) 
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where 

B{(p)n{p',j',k')n{p,j,k) = Sppf exp{i2']T{N'^jj' + M'^kk')} (3.3) 

^(0)n(l,jl,fci)n(2j2,fc2) = ~M'P)n{2,J2,k2)n{l,ji,ki) = 

tanh(3Eiexp{t2TTN-'ji}6{ji,j2,k,,k2) (3.4) 

-4(0)„(3ji,fci)n(4,J2,fc2) = ~M(P)n(4,J2,k2)n(3,ji,ki) = 

tanh (3E2exp{i27rM-^ki}6{ji, j2, ki, k2) (3.5) 

^{4>)n{l,ji,ki)n{3,J2M) = ^^i4>)n(3,J2,k2)n(l,ji,ki) = 
^(^)'i(lji,fei)n(6j2,fc2) = -M(P)n(6,J2,k2)n(l,jiM) = 
-^(0)n(2,ji,fci)n(4,J2,A:2) = ^(0)n(4,J2,fc2)n(2,ji,fci) = 
^(0)n(2,ii,A:i)n(5,i2,A:2) = " ^(0)n(5,J2,fc2)n(2,ji,fci) = 
-M4>)ni3Ji,kiM6,J2,k2) = ^(0)n(6,J2,fc2)n(3,ii,fci) = 
-^(0)n(4,ji,fci)n(5,J2,fc2) = ^(0)n(5,J2,fc2)n(4,ii,fci) = 
^(0)n(5Ji,A;i)n(6j2,A;2) = ^^(0)n(6,J2,fc2)n(5ji,fci) = 

5iJi,J2,ki,k2) (3.6) 

SiJij2, ki, k2) = {MN)~-^{5N-n,J2 + ^hN5j2N){hi-k^,k2 + h^M5k2M) (3.7) 

All the other elements of the matrix A(0) are equal to zero. 
Due to Proposition 1 from ([|12], chap. IX, sect. 5) 



Pf {B^A{(t))B) = det 5Pf i(0). (3.8) 



In view of the definition ( p.3|) the matrix B = I^® Cn ® Cm, where Iq is the identity 
6 X 6 - matrix and N x N - matrix {CN)kj = exp{i27rN~^kj}. Therefore 

det B = (det Civ)^*^(det CmT^- (3.9) 

For the determinant of the matrix {CN)jk = exp{i27TN~^{N — j)k} we have 

det Cn = (-1)'^' det Cn (3.10) 

where [r] denotes the integral part of the real number r. The relation (|3.1|) implies 

CnCn = NIn (3.11) 

where I^ is the identity N x N - matrix. It follows from the relations (|3.1(]| ), ( p.ll| ) 



rJV-li 



(detC^)2 = (-l)[^lAr^. (3.12) 

The substitution of the equality (|3.12|) into the equality (|3.9| ) yeilds 



detS = {-iyn^i+^'^i^\MNy'''\ (3.13) 
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Instead of the definition ( |2.8| ) we make use of the definition from ([0, chap. IX, sect. 
5). Let Cfc, k = 1, ..., 6MN, be the basis of the space where the matrix A(0) acts. Then 

(6MN _ \ 

Y^ i(0),fce, A e J = 23*^^(3MiV)!Pf i(0)ei A ■ ■ ■ A e^MN. (3.14) 

j,k=l J 

We introduce 6x6- matrix 

{D{j,k))i,i2 = M N A{(t))n(i^,j,k),n(i2,N-j,M-k) 
{D{N,k))i^i^ = MNA{(f))n(i^^N,k),n{i2,N,M-k) 
{D{j,M))i^i2 = MN A{(f))n(i^j^M),n(i2,N-j,M) 

{D{N, M))i,i, = MNA{<p)n(,^^N^M),n(i2,N,M) (3-15) 

where 1 < «i < 6, 1 < ^2 < 6, 1 < j < A^, 1 < A; < M . In view of the definitions (|3T 



the matrix (|3.15|) is anti - Hermitian: {D{j,k))i^i-^ = —{D{j,k))i-^i^. For the natural 



numbers j = y, iV, k = ^^^M the matrix -D(j, k) is real and therefore antisymmetric. Hence 
its Pfaffian is defined for these natural numbers j = ^^N ^ k = ^^^M 

FiD{j,k) = 
Z1Z2 exp{i2n{N-^j + M'^k)} + zi exp{i27riV"^j} + Z2 exp{i27rM-^A;} - 1 (3.16) 

where Zi = ta.nh.f3Ei, i = 1,2. The determinant of the matrix D{j,k) is defined for any 
numbers 1 <j < N,l <k< M 

det D{j,k) = 
(1 + zl){l + zl) - 2zi{l - z^) cos{2ttN-^j) - 2^2(1 - zf) cos{27rM-^k). (3.17) 



The equality ( |1.4| ) coincides with the relation 

det D{N, M) = {FiD{N, M)f . (3.18) 



By making use of the definitions ( p.4| ) - ( |3.7| ) and the Pfaffian definition ( p. 141 ) it is possible 



to show that the Pfaffian of the QMN x 6MN - matrix A((j)) is proportional to the product 
of all Pfaffians of 6 x 6 - matrices D{j, k) for the natural numbers j = ^,N , k = ^^,M and 
the product of all determinants of the matrices D[j,k), I < j < N, 1 < k < M, except 
j = ^,N , k = ^-,M simultaneously. 
Let us introduce the products 

I2 rn2 

Pu{h,l2;m,,m2)= n n [i^ + 4)i^ + 4)- 

j=li k=m\ 

-2zi{l - zl) cos{2nN~^j) - 2^2(1 - ^i) cos(27rM"^A;)] (3.19) 

h rn2 

P2i{h,l2;m,,m2)= n n [(l + -2i)(l + ^2)- 

j=h k=mi 

-2zi{l - zl) cos{nN-\2j - 1)) - 2^2(1 - zf) cos{27r M-'^k)] (3.20) 
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Pi2(/i,/2;mi,m2)= n n [(l + -2i)(l + ^2)- 

j=h k=mi 

-2zi{l - zj) cos{27rN-^j) - 2;Z2(1 - zf) cos{7rM~\2k - 1))] (3.21^ 



P22(/i,/2;mi,m2) = n n [(i+^?)(i+-22)- 

j=ll k=m\ 

-2zi{l - zl) cos{7tN-\2j - 1)) - 2^2(1 - zl) cos{nM-\2k - 1))]. 

Calculating the Pfaffian Pf A(0) and making use of the relations ( |3.2| ), ( |3.8|) and 
we have: 
if the numbers A^, M are odd, then 

Pf A(0) = Fn{zi, Z2){z^Z2 + Z1 + Z2-I), 

where Zi = tanh/5i?j, i = 1,2 and 

, N -1 M -1, , N -1 M-\, 
Fnizu Z2) = Pii(0, ^— ; 1, -^—)P,,(1, -^-; 0, -^—); 

if the number A^ is odd and the number M is even, then 

Pf v4(0) = ^12(2:1, Z2){-~ZiZ2 + Zi- Z2- l){ZiZ2 + Zi + Z2 - 1), 

A^-1 M N -I M 

F,2{zu Z2) = Pn(0, ^-; 1, y - l)Pn(l, ^-; 0, -); 

if the number A^ is even and the number M is odd, then 

Pf A(0) = ^13(2:1, Z2){-ZiZ2 - Zi+ Z2- l){ziZ2 + Zi + Z2 - 1), 

A^ M-1 A^ M-1 

Fi3(^i, ^2) = Pn(0, -; 1, ^— )i^ii(l, y - 1; 0, ^— ); 

if the numbers A^, M are even, then 

Pf A(0) = Fi4(zi, 2;2)(2;i2;2 - zi - ^2 - l)(-2i22 - 2;i + 2:2 - 1) X 

{-Z1Z2 + Zl- Z2~ l){ZiZ2 + Z1 + Z2- 1), 

N M N M 

FiM, Z2) = Pii(o, y ; 1, Y - i)Ai(i, Y - 1' 0' y)- 



Analogously the definitions ( f^.49| ) - ( |2.57| ) imply: 



if the numbers A^, M are odd, then 

Pf A(0 + a*) = F2i{Zu Z2){-ZiZ2 -Z1+Z2-I), 

if the number A^ is odd and the number M is even, then 

Pf y4(0 + a*) = F22{Zi, Z2){ziZ2 - Zi - Z2 - l){-ZiZ2 - Zi + Z2 - I] 
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3.22) 



3.23) 

3.24) 

3.25) 
3.26) 

3.27) 
3.28) 



3.29) 
3.30) 



3.31) 
3.32) 

3.33) 



A^ + 1 M N -I M 

F,,{z„ z,) = P2i(l, ^; 1, y - l)P2i(l, ^-; 0, - 

if the number N is even and the number M is odd, then 

, A^ M-1 , A^ M-1 
P2i(l, y ; 1, ^— )P2i(l, y ; 0, ^— ); 

if the numbers A^, M are even, then 

PfA(0 + a*)=F24(^l,^2) = 

N M N M 

P2i(l,y;l,y)P2i(l,y;0,y-l). 



Similarly the definitions (3) - (jIIS^ ), (|3|) - (|2:56|) and (|38D imply: 
if the numbers iV, M are odd, then 

Pf A(0 + b*) = P3i(2;i, Z2)(-2l2;2 + 2;i - ^2 - 1), 

, N-1 M-1 , iV-1 M + 1, 
^31(^1, ^2) = Pi2(0, ^— ; 1, ^— )Pi2(l, ^— ; 1, ^— ); 

if the number A^ is odd and the number M is even, then 

PfA(0 + 6*)=P32(^l,^2) = 

, N-1 M,^ , N ~\ M, 
Pi2(0, -^—, 1, y)Pi2(l, -^—, 1, y ); 

if the number A^ is even and the number M is odd, then 

FiA{(p + b*) = F^si^zi, Z2){ziZ2 - Zi - Z2 - l){-ZiZ2 + zi - Z2 - 1), 

A^ M-1 A^ M+1 

P33(^i, Z2) = Pi2(0, y ; 1, ^— )^i2(l, y - 1; 1, ^p); 

if the numbers A^, M are even, then 



PfA(0 + 6*)=P34(^l,^2) = 



PuiO 



N 



M. 



N M, 



2 i;i'y)^i2(i,y;i, 2 



Similarly the definitions (|^ - (jIIS^ ), (|3|) - d^M) . (HI), (HO) imply: 
if the numbers A^, M are odd, then 

FiA{(f) + a* + b*) = Fii{zi, Z2){ziZ2 - zi - Z2 - 1), 
if the number A^ is odd and the number M is even, then 

PfA(0 + a* + 6*)=P42(^l,^2) = 

N + 1 M N-1 M 

"22(1, — ^ — ; 1, y)-P22(i, — ^ — ; 1, y); 



(3.34) 



(3.35) 



(3.36) 



(3.37) 
(3.38) 



(3.39) 

(3.40) 
(3.41) 



(3.42) 



(3.43) 
(3.44) 



(3.45) 
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if the number A^ is even and the number M is odd, then 



A^ M-l N M+1 

"22(1, —] 1, — ^ — )-P22(l, —] 1, — ^ — ); 



(3.46) 



if the numbers A^, M are even, then 



ViA{(t) + a* + h*)=F^i{z^,Z2) = 
(P22(l,y;l,y)f. 



(3.47) 



The formulae ( 3.23| ) - ( 3.47| ) are in accordance with the formulae for detAj. The 
formulae i^, (lOSi) - (|^ imply 



detA(0) = Pn(l,Ar;l,M) 

det A{<p + a*) = P2i(l, A^; 1, M) 

detA{4> + b*) = Pi2{l,N;l,M) 

det A(0 + a* + b*) = P22(l, A^; 1, M). 



The expressions (|3.48| ) coicide with the expressions for dety44, i 



(3.48) 
, 4, given in the book 



4 Free Energy 



The following lemma is the staightforward consequence of the relation (|1.4|) . 
Lemma 4.1 For the variables Zi = tanh j3Ei, j3 > 0, Ei ^ 0, i = 1,2 the polynomial 

(1 + zl){l + zl) - 2Zi{l - zl) COS01 - 2^2(1 - z\) COS02 > 

ij at least one of the followng conditions 

cos (pi = sgn Zi,i = 1,2, 

is not valid. 

Proof. Since < \zi\ < 1, i = 1,2, the following inequality holds 

{l + zl){l + zl) - 2^1(1 -Z^) cos 01 -2^2(1 -^i)cOS02 > 

(1 + zf){l + zl) - 2|^i|(l - zl) - 2\Z2\{1 - zl) 



Due to the equality ( |1.4| ) the right hand side of the inequality 

{\ZiZ2\ + l^ll + \Z2\ - if. 



is equal to 



(4.1) 
(4.2) 



(4.3) 



(4.4) 



The inequality ( [4.3|) may be an equality if and only if both conditions (|4.2| ) are valid. The 
lemma is proved. 

Lemma 4.2 The range of values for the variables Zi = tanh f3Ei, f3 > 0, Ei ^ 0, i = 1,2 is 
divided into eight two dimensional domains: 



[1 + zi){l + Z2) > 2, (1 - zi){l + Z2) < 2, 

'1+Zi){l-Z2)<2,{l-Z^){l-Z2)<2, 

Zl > 0, Z2 > 0; 



(4.5) 
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(1 + ^1 
;i + ^i 



[1 + ^1 

'l+Zi 



[l + Zi 
'l + Zi 



[l + Zi 
'l + Zi 



Z\ 



[1 + ^1 
;i + ^i 



(1 + ^1 
;i + ^i 



and into four one dimensional 



Zl 
Z\ 



[1 + ^1 
;i + ^i 



[1 + ^1 
;i + ^i 



)(l + ^2) 


<2, 




-^i)( 


)(l-^2) 


<2, 




^1 


)(l + ^2) 


<2, 




-^i)( 


)(l-^2) 


<2, 




-^i)( 

-^1 


)(l + ^2) 


<2, 




-^i)( 


)(l-^2) 


<2, 






)(l + ^2) 


<2, 




-^l)( 


)(l-^2) 


>2, 




-2l 


)(l + -22) 


<2, 




-^l)( 


)(l-^2) 


<2, 




Z\ 


)(l + ^2) 


<2, 




-^i)( 


)(l-^2) 


<2, 




•^1 


)(l + ^2) 


<2, 




-^i)( 


)(l-^2) 


<2, 




-2l 


/ domains: 






)(l + -22) 


= 2, 




-^l)( 


)(l-^2) 


<2, 




Zl 


)(l + ^2) 


<2, 




-^l)( 


)(l--22) 


<2, 




-^l)( 

-^1 


)(l + ^2) 


<2, 




-^l)( 


)(l-^2) 


= 2, 




-^l)( 



1 + ^2) < 2, 

1 - ^2) < 2, 

> 0,^2 > 0; 



- + -22) > 2, 

- - ^2) < 2, 
< 0,^2 > 0; 



1 + ^2) < 2, 
1 - ^2) < 2, 

< 0,^2 > 0; 

1 + ^2) < 2, 
1 - ^2) < 2, 

> 0,^2 < 0; 

1 + ^2) < 2, 
1 - ^2) < 2, 

> 0,2:2 < 0; 

1 + ^2) < 2, 
1 - ^2) > 2, 

< 0,^2 < 0; 

1 + ^2) < 2, 
1 - ^2) < 2, 

< 0,^2 < 0; 



1 + ^2) < 2, 

1 - ^2) < 2, 

> 0,^2 > 0; 

1 + ^2) = 2, 

1 - ^2) < 2, 

< 0,2:2 > 0; 

1 + ^2) < 2, 
1 - ^2) < 2, 
2:1 > 0, 2:2 < 0; 



(4.6) 



(4.7) 



(4i 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 
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1 + z{){l + Z2) < 2, (1 - z{){l + Z2) < 2, 
l + Zi)(l-^2)<2,(l-^i)(l-^2) = 2, 

2i<0,Z2<0. (4.16) 

1,2 be given. Let us prove that if 



Proof. Let two sets of numbers ej = ±1, ctj = ±1, z = 
these two sets do not coinside, then two inequahties 

(l + eiZi)(l + e2Z2)>2 
[l + aizi){l + a2Z2) > 2 



are not compatible for the variables z, 
eiCTi = 



tanh/3i?j, i 



(4.17) 
1,2. Indeed, let ei 7^ ai, i. e. 



— 1. Multiplying the inequalities ( |4.17| ) we obtain the inequahty 

[l - Zl){l + e2Z2){l + <J2Z2) > ^ 

Since zi = tanhpEi, (3 > 0, Ei ^ 0, i = 1,2, then < \zi\ < 1 and 

0<1-Zi<l,0<l + e2Z2 < 2, < 1 + a2Z2 < 2. 



(4.18) 



(4.19) 



The inequalities ( [4.1(j| ) and ( |4.19| ) are not compartible. The case 62 7^ ^2 is considered 
similarly. 

Let us prove that for the variables Zi = tanh j3Ei, (3 > 0, Ei y^ 0, i = 1,2, the first 
inequality (f4.17| ) implies eiZi > 0, i = 1,2. Indeed, for the variables Zi the inequality 
< 1 + eiZi < 2 is valid. Hence the first inequality ( [4.17|) implies 1 + eiZi > 1 or ejZj > 0. 
The lemma is proved. 
Theorem 4.3 For the partition function ( p.2|) the following relation is valid 



p-^ lim {MN)-HnZ{N,M) 

M,N^oo 



r2TT 



-/5-Mn(2 cosh /?Ei cosh /3E2) - l/2/3-^(27r)-^ / "rf^i / (102 

Jo Jo 

ln[(l + zl){l + zl) - 2zi{l - zl) cos 01 - 2^2(1 - zl) cos02] 



(4.20) 



where the variables Zi = tanh (3Ei, (3 > 0, Ei y^ 0, i = 1, 2. 

Proof. Let the numbers A^, M be odd. Then the homological formula ( |2.42| ) and the equalities 

O), (pTD, (p3), (HD imply 

Z{N,M) = {2 cosh (3 El cosh f3E2)^^ x 
l/2[Fn(^i, ^2)((1 + ^i)(l + Z2) - 2) - F2i(zi, Z2)((l - ^i)(l + Z2) - 2) - 

-F3i(2i, ^2)((1 + ^i)(l - Z2) - 2) - F,i{zi, ;22)((1 - 2i)(l - Z2) - 2)]. (4.21) 

Since the function cosx is monotone decreasing on the interval [0, vr], we have 

A^-1 



cos(27riV-^j) < cos(7riV-^(2j - 1)), 1 < j < 
cos(7rA^~^(2(j + 1) - 1)) < cos(27riV"^j),0 < j < 



2 ' 

iV-3 



(4.22) 



It follows from the relations ( ^:24D , (g^D, ( ^381) , {^A§ and the inequalities (j^) that for 

Zl > 0, Z2 > 



N+1 N+1 

-r2i(, ;; , , i, 



^21(^1,2^2) < 
^— )(Pn(0, 0; 1, ^—))-^Fiiizi, Z2) 



(4.23) 
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f,.„„,,) > P,.,U;0,^)P.(^^.^^;l,^i^i X 



2 ' 2 ' ' 2 



,„ ,N-l N-l M-l,„ ,A^-1 iV-1 M-L, , , 



2 ' 2 



2 ' 2 



(4.24) 



P.,(l, ^; £ + i, ^)(P„(1, ^; 0, Oj-Fl'd I) 



(4.25) 



^3.,,, ,) > P.,(0, i^; , 1)P.,,1, ^; ^, ^) . 

, , A^ - 1 Af - 1 M - 1 , iV - 1 M - 1 M - 1 , 1 

(Pn(0, ^— ; ^— , ^— )Pn(l, ^— ; ^— , ^— ))-iFn(^i, ^2) (4.26) 



^ , , „ ,, AT-l M+1 M+l,„ ,iV+l AT + l ^ M-1. 
^^41(2^1,22) <P22(1, ^^; „ , ^ )P22(^^.^^;1. — z — ) X 



' 2 ' 2 ' 2 



2 ' 2 



(Pn(l, ^^^; 0, 0)Pn(0, 0; 1, ^^LA))-^ F^^{z^, z^) (4.27) 



A^- 1 M- 1 

P4i(^i, ^2) > ^22(1, ^— ; 1, 1)^22(1, 1; 1, ^— ) X 

^ ,^ N-l M+l M + 1,^ .A^+l A^ + l ^ M-1, 
i^22(l, ^— ; ^— , ^— )P22(^— , ^-; 1, ^— ) X 

, , iV-3 M-1 Af-1 ,iV-l iV-1 M-3,, 1 
(Pii(0, ^— ; ^— , ^— )Pn(^— , ^— ; 0, ^— ))-' X 

,^ , A^-1 M-1 M-r^ ,A^-1 A^-1 M-L, , , , , , 

(Pn(l, ^— ; ^— , ^— )Pn(^— , ^— ; 1, ^— ))-^Pn(;^i, z^). (4.28) 



The relation (^4.21 ) and the inequahties ( 4.23|) , ( [4.25| ), (^4.27 ) imply for the domains 
D and (|41^ ) 



2Z(Ar,M)(2cosh/5Eicosh/5E2)-^^^(Pii(zi,22))"^ < 

((l + Zi)(l + Z2)-2)- 

-((1 - .0(1 + ~~.) - 2)/^.(^. ^; 1. ^)(Pn(0. 0; 1. ^))- - 
-((1 + .)(! - ..) - 2)P.(1. ^; ^, ^)(P„(1, ^; 0, 0))- - 

-((1 - 2;i)(l - Z2) - 2)P22(1, ^— ; ^— , ^— )^22(^— , ^— ; 1, ^— ) X 

(Pn(l, ^^; 0, 0)Pn(0, 0; 1, i!l^))-i(4.29) 



By making use of the inequalities ( |4.22|) we have 



A^- 1 
Pn(l,^^;0,0)(P22(l, 



iV- 1 M+l M + l ,_i 

)) < 



2 ' ' /^ --^ ' 2 2 2 

((1 + zl){\ + ^2) + 2^i(l - 4) cos(7riV-^) - 2^2(1 - z\)) x 
((1 + z\)(\ + z^) - 2^1(1 - zl) co^(txN-^) + 2z2(l - ^i))"' 



(4.30) 
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Ai(0,0;l. 



^'-l),p,,(£±i,^;l,^))-.< 



2 /v -V 2 ' 2 2 

((1 + zl){l + zl) - 2^1(1 - zj) + 222(1 - zf) cos(7rM-i)) x 

((1 + zf){l + zl) + 2^1(1 - zl) - 2z2(l - z\) cos{nM-')y^ 



(4.31) 



iV + l iV + 1 M-1 iV + l iV + 1 M-1 

-^2ll ^ , ^ , i, ^ Jl-r^22l ^ , ^ ,-L, ^ JJ < 

((1 + zf){l + Z2') + 2zi{l - zl) + 2^2(1 - zl) cos{ttM-^)) X 
((1 + zf){l + zl) + 2zi{l - zl) - 2^2(1 - zf) cos{7rM-^))-' 



(4.32) 



iV-l_M + l M + 1 iV-l,M + l M+1 

-r^l2U, — ^ — , ^ , ^ Jl-r22U, — ^ — , ^ , ^ jj < 

((1 + zl){l + zl) + 2z,{l - zl) cos(7rAr-^) + 2^2(1 - %')) x 
((1 + zf){l + zl) - 2^i(l - zl) cos{7tN'^) + 2^2(1 - ^i))"'- 



(4.33) 



In view of Lemma 4.1 the right hand sides of the equahties (|4.30| ) - ([4.33|) are bounded. 
By dividing the inequailty ( |4.29| ) by 

^ ,^ N-1 M+1 M+l.^ .N + 1 N + l ^ M-1, 
^22(1, ^7^—; — ^ — , —^ — )i'22(^7^, ^^; 1, ^x 



(^11(1 



2 ' 2 

A^-l 



2 ' 2 
;0,0)Pn(0,0;l, 



2 '-"/'-v--- 2 
and by using the inequahties ( [4.30| ) - ( |4.33|) we have 



"-'rr^ 



-p-^ hm {MN)-HnZ{N,M)> 

M.N^oo 



-f3-Hn{2cosh(3Eicoshf3E2) - l/2/3-^(27r)-2 / "rf^i / " rf^2 

JO Jo 

ln[(l + zl){l + zl) - 2zi{l - zl) cos^i - 2^2(1 - zf) cose2] 



(4.34) 



for the domains ([4.5|) , (4.6) and ( 4.13 ). 

The relation ( [4.21| ) and the inequahties ( [4.24| ), ( [4.26| ), ( [4.28| ) imply for the domains (^ ) 
(PD and (|41^ ) 



2Z(Ar,M)(2cosh/?Ei cosh /3E2)-^-^^(Fn (21,22))"^ > 

iil + z,){l + Z2)-2)- 



-((l-Zi)(l + ^2)-2)P2i(l,l;0 



M-1 



lP2l( 



A^ + 1 iV + 1 ^ M-1 



:1. 



(P: 



,N-1 iV-1 „ M-1 



111 



-((i + 2i)(i-22)-2)Pi2(a 



iV-l 



iP 
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. iV-l iV-1 .^ M-l ^^_i 

^ 2 ' 2 



' ' 2 



-))- 



;i,i)i^i2(i, 



A^-l M+1 M + 1, 



2 ' 2 2 ' 

A^-1 M-1 M-1 A^-1 M-1 M-1 ^1 

("11(0, — - — ; — - — , — - — )Pii(l, — - — ; — - — , — - — )) 

A^- 1 M-1 

-((1 - ^0(1 - ^2) - 2)P22(1, ^— ; 1, 1)P22(1, 1; 1, ^r ' 



X 
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N-1 M + 1 M+1 N + 1 N + 1 M-1 

-r22i-L, — - — ; — - — , — - — j-r22i — - — , — - — ; i, — - — ) X 
(Pn(0. ^; ^, ^,P„(^. ^; 0. ^))- X 

(Pnd. ^■. ^, ^)Pn(^. ^. 1. V))-'. (4.35) 



By using the inequalities (|4.22|) it is easy to obtain for the domains ( [4 .51) , (|]6|), (|4.13| ) 
for large numbers A^, M the inequalities similar to the inequalities ( [4.3CI| ) - ( |4.33| ) for the 
products 



[Pu 



P2i(i, 1; 0, ^— )P2i(^— , ^— ; 1, ^— )x 

,A^-l A^-1 M-1 „ ,N -\ N -\ M-1 



2 ' 2 



-;0, 



)Pi 



))-\ 



, iV-1 ,„ , iV-1 M+1 M+1, 
Pi2(0, ^— ; 1, l)Pi2(l, ^— ; ^— , ^^)x 



, , N-\ M-1 M-1 , iV-1 M-1 M-1 , 1 
(Pn(0, ^— ; ^— , ^— )Pn(l, ^— ; ^— , ^— ))-\ 



2 ' 2 ' 2 

AT-l 



2 ' 2 



Po.fl, 



M-1. 
;1,1)P22(1,1;1,— )x 



P..fl 



2 ' ' ^ --^ ' ' ' 2 
A^-1 M+1 M+ 1 „ ,iV + l iV + 1 M-1, 



221J-, 



(Pii(0, 

(Al(l 



)P22(- 

)Pxx\ 



2 ' 2 ' 2 ''""' 2 ' 2 ' ' 2 
iV-3 M-1 M-l.„ .N -\ N -\ n^~^ 



2 ' 2 ' ' 2 



2 ' 2 ' 2 

A^-1 M-1 M-L„ .iV-1 iV-1 . M-1 



))-ix 



)Pii( 



2 ' 2 ' 2 '''"' 2 ' 2 ' ' 2 
Therefore the inequality ( [4.35| ) implies the following inequality 



)V. 



p-^ lim {MN)~HnZ{N,M) < 

M,N~*oo 



27r f2Tr 



-(3-Hn{2cosh(3Eicosh(3E2) - l/2l3-\2n)-^ f ^ d0i f 

Jo Jo 



d9o 



ln[(l + zf){l + z^) - 2zi{l - zl) cos^i - 2^2(1 - zl) cos02 



(4.36) 
(EH) 



for the domains (|]5D, (U), (|]T|). The equality (|]20|) for the domains (|]5|), 
follows from the inequalities (|4.34| ), ( [4.36| ). 

Let the numbers N,M be odd and z^ < 0, ^2 > 0. The relations (|;2§, ( ^32]) , (|338|) 
( p.44|) and inequalities ( [4.22|) imply 



i^ii (^1,2:2) < 

M-1 , ,A^ + 1 A^+1 M-1 , 1 
Pn(0, 0; 1, ^_)(P2i(^— , ^— ; 1, _^))-iF2i(zi, ^2) 



(4.37) 
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Fu{Zi,Z2) > Pu 



.N-1 N-1 M-1. 



2 2 2 ' 

,N+1 N + 1 M-1 

V 7^ ' 7^ ' -'-5 7^ 



Pn 



.N-1 N-1 M-1. 



{P2i{'-^r^, '-^■, 1, ^^— ^)P2i(l, 1; 0, 



2 2 

M-1 



))-'F2iiz,,Z2) (4.38) 



M-r„ , A^-1 M + 1 M+r 
F3i(^i, Z2) < Pi2(0, 0; 1, ^— )Pi2(l, ^— ; -^—, -^—) X 

. . N-1 , ,iV+l A^+1 M-1 , -I 

(P2i(l, ^^; 0, 0)P2i(^^, ^^; 1, ^— ))-^F2i(;2i, ^2) 



(4.39) 



. A^-1 ,^ , iV-1 M+1 M+1, 
F3i(;^i,^2)>P2i(l,^— ;0,0)Pi2(l,^— ;^— ,^— )x 

, N-1 N-1 M-1 ,., , N-1 M-1 M-1,, . 

(Pi2(^^, -^^■, 1, ^— ))'(P2i(i, -^^■, -^^, -^^)r' X 



fP2i(l,l;0,^^^)P2i( 



2 ' 2 

1, 



N+1 N + 1 . M-1 1 

;^ — )) F2i{zi,Z2) 



(4.40) 



, A^-1 M+1 M + 1 
^41(^1, Z2) < ^22(1, ^— ; ^— > ^— ) ^ 

(P2i(l,^^;0,0))-iP2i(^i,^2) 



(4.41) 



, iV + 1 , , N-1 M+1 M+1, 

P4l(^l, Z2) > P2l(l, ^— ; 0, 0)P22(1, -^—, -^—, -^—) X 

, , N+1 M-1 M-1, , N -1 M -1 M -1,, , 

(P2i(l, ^— ; ^— , ^— )P2i(l, ^— ; ^— , ^— ))-^P2i(^i, Z2). (4.42) 

By making use of the relation ( [4.21| ) and the inequahties ( |4.38|) , ( |4.39| ), ( [4.41|) it is possible 
to obtain the inequality ( [4.34|) for the domains ( [4.7|) , ([4.8|) , (|4.14| ). Similarly the relation 
( ^^211) and the inequalities ( ^^ , (|^, (Fl42| ) imply the inequality (14361 ) for the domains 
0, (U), (|4A^ ). The inequalities (|3|), ( CT) for the domains Q, (|t|), (|A|) imply 
the equality ( |4.2(J| ) for these domains. 

The proof of the equality ( |4.2(]| ) for odd numbers N, M and for the domains ( [1.9|) , ( ^l.lOp , 
( [4. 151 ) is quite similar to the proof of the equality for odd numbers N, M and for the domains 

(O), (pp, (gj§- 

Let the numbers N,M be odd and zi < 0, Z2 < 0. The relations (|]2|), ( ^321) , (^M) , 
( p.44| ) and the inequalities ([4.22|) imply 



A^- 1 M-1 

Pn(^i,^2)<Pii(l,^— ;0,0)Pn(0,0;l,^— )x 

,P,,(1. E^, £±i, ^,P,,(^, ^; 1, ^r^F^M,.) 



Fu{zi,Z2)>Puil, 



2 2 ' ' 2 

A^-1 M-1 M-1, 



(4.43) 



Pii 



2 2 2 ■ 

iV-1 iV-1 M-1 



„ ,„ Ar-3 M-1 M-1, 
Pii(0,— ^;^ — ,^ — ) X 



-;i, 



Pii 



2 2 2 ' 

,A^-l A^-1 M-3. 



2 ' 2 

(P22(1,^^;1,1)P22(1, 

,A^+1 A^+1 M- 1, 



^ 2 ' 2 "' ' 2 ''^ 

A^-1 M + 1 M+ L,_i 

)) X 



{P22{'-Tr^, '-^■, 1, ^^^,— )P22(1, 1; 1, 



2 ' 2 
M-l._i 



))-^P4i(^i,^2) (4.44) 



23 



, N -\ ,, , N -\ M + 1 M + 1 , 1 , , , , 

P2l(l, ^— ; 0, 0)(P22(1, ^— ; ^— , ^— ))-^F4l(Zl, Z2) (4.45) 

, iV-lM-lM-1 , A^+lM-lM-1, 

i^2i(^i,^2) >P2i(l,^— ;^— ,^— )Ai(i,^— ;^— ,^— ) X 

(P,,,, ^; 1, l)p,,(l, ^; £±i, ^,)-.p„(,, .., ,4.46, 

P.(0, 0; 1, ^)(P.,i^, ^; 1. ^)rC(l .!) (4.47) 

,N ~1 N -\ M-1 ,Ar-l AT-l M + 1, 
^31(2^1, Z2) > Pi2(^— , -^—, 1, ^— )Pi2(^— , ^— ; 1, -^—) X 

(P22(l, 1; 1, ^— )P22(^-, ^— ; 1, ^— ))-^P4i(^i, ^2). (4.48) 

The inequalities ( |4.43| ) - ( |4.48| ) imply the inequalities ( [4.34| ), ( [4.36| ) and therefore the equality 
( ^:20D for the domains (|;T1|), ( ^l2D , ( ^+61) . The equality (|O0|) is proved for the odd 
numbers N, M. 

The cases when N is odd, M is even; when N is even, M is odd and when A^, M are even 
are analogous to the case when the numbres N, M are odd. The theorem is proved. 
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